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Abstract 

Few-body problems involving Coulomb or gravitational interactions between pairs of particles, 
whether in classical or quantum physics, are generally handled through a standard multipole ex- 
pansion of the two-body potentials. We develop an alternative based on an old, but hitherto 
forgotten, expression for the inverse distance between two points that builds on azimuthal sym- 
metry. This alternative should have wide applicability throughout physics and astronomy, both 
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for computation and for the insights it provides through its emphasis on different symmetries and 
structures than are familiar from the standard treatment. We compare and contrast the two meth- 



ods, develop new addition theorems for Legendre functions of the second kind, and a number of 
useful analytical expressions for these functions. Two-electron "direct" and "exchange" integrals in 
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many-electron quantum systems are evaluated to illustrate the procedure which is more compact 
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than the standard one using Wigner coefficients and Slater integrals. 
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I. INTRODUCTION 



For pairwise Coulomb or gravitational potentials, one often expands the inverse distance 
between two points x and x' in the standard multipole form [[I ] 
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where r< (r>) is the smaller (larger) of the spherical distances r and r', and P^(cos7) is the 
Legendre polynomial with argument 

cos 7 = x ■ x' = cos 9 cos 9' + sin 9 sin 9' cos(0 — </>')• (2) 

The set of six coordinates {x, x'} may be viewed either as defining two points relative 
to the origin or as the coordinates of a three-body system once the motion of the center 
of mass has been separated. In the "body frame", three out of the six coordinates are 
dynamical variables, the potential energy depending only on them. Of various choices for 
these variables, one is the set of three separation distances, another the triad (r<, r>, 7) as in 
Eq. ([[]). With respect to a space-fixed "laboratory frame," three more angles constitute the 
full set of six coordinates, the choice in Eq. (H) of (9, 9', — 0') being suited to the spherical 
polar coordinates of the individual vectors; thus, x: (rsin#cos0, rsin#sin0, rcos9). 

In this paper, we present an alternative expansion to Eq. (P based on cylindrical (az- 
imuthal) symmetry which should be of wide interest in physics and astrophysics. This ex- 
pansion arose in a recent investigation by two of us [[|] of gravitational potentials in circular 
cylindrical coordinates x = (R, 0, z): 
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with Q m _i a Legendre function of the second kind of half- integer degree H, and x defined 
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as 



_ R 2 + R' 2 + - z' f _ r 2 + r' 2 - 2rr' cos 9 cos 9' 
X= 2RR' ~ 2rr'sin^sin 9' ' ^ ' 

Although this expansion has been recorded in many places |5], § , its full significance has 

not been appreciated. We have traced its earliest occurrence to the work of E. Heine in 

the mid-nineteenth century || and will, therefore, call it the "Heine identity." At its most 

general, it takes the form 
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reducing to Eq. ([3D when applied to the inverse distance between two points. We have found 
it to be a more efficient approach for problems with a cylindrical geometry and have used it 
for compact numerical evaluation of gravitational potential fields of several axisymmetric and 
nonaxisymmetric mass distributions We now set Eq. (|3[) in a broader context, together 
with new associated addition theorems and a novel application in quantum physics, hoping 
to encourage wider use of this expansion throughout physics. 



II. THE ALTERNATIVE EXPANSIONS 

The expansion in Eq. (p]) disentangles the dynamics contained in the radial variables 
from symmetries, particularly under rotations and reflections, pertaining to the angle 7. 
Whereas the three variables {r<,r>,7} at this stage are joint coordinates of x and x', a 
further disentangling in terms of the independent coordinates so as to handle permutational 
and rotational symmetry aspects of the problem is often useful and achieved through the 
addition theorem for spherical harmonics 
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where Yg m are the usual spherical harmonics Jj]], T(z) is the gamma function, and P^i^z) is 
the integer-order, integer-degree, associated Legendre function of the first kind [[|. Using this 
to replace P^(cos7) in Eq. ([]]), we obtain the familiar Green's function multipole expansion 
in terms of all six spherical polar coordinate variables x and x', 



|x-x'| ^f^K^J m tl^r(£ + m + l) 
Apart from the first factor with dimension inverse-distance formed from the geomet- 
ric mean of the two lengths r and r', this expression involves only four combinations: 
(r</r>, 9, 8', <ft — <$') of the six coordinates x and x'. In spite of widespread familiarity with 
the multipole expansion, this reduction from six to four essential variables has not been 
appreciated fully. We were led to it by the parallel investigation below of the alternative 
expansion, and note that it is a natural consequence of the separation distance being inde- 
pendent of the orientation of that separation in the laboratory frame, and thus independent 
of two angles serving to specify that orientation. 



This multipole expansion in terms of spherical harmonics is very broadly utilized across 
the physical sciences. For example, with i and m interpreted as the quantum numbers of 
orbital angular momentum and its azimuthal projection, respectively, a whole technology 
of Racah-Wigner or Clebsch-Gordan algebra is available || for handling all angular (that 
is, geometrical or symmetry) aspects of an iV-body problem, the dynamics being confined 
to radial matrix elements of the coefficients (r < /r > ) £+ ^ in Eq. ([I]). Although many other 
systems of coordinates have been studied for problems with an underlying symmetry that is 
different from the spherical, Eq. (|l|), in combination with the addition theorem for spherical 
harmonics, has gained such prominence as to have become the Green's function expansion 
of choice even for non-spherically symmetric situations. 

But, consider now the expansion in Eq. (^), which may be viewed either, in analogy with 
Eq. (|IJ), as an expansion in Legendre functions, now of the second kind in the joint variable 
X of the whole system, with coefficients {RR')~ l / 2 e im ^~^'\ or as a Fourier expansion in the 
variable (0 — <$') with the Q's as coefficients. In this latter view, a further step allows us to 
develop a new addition theorem for these Legendre functions. Interchanging the i and m 
summations (Fig. in Eq. (0), we obtain 

1 S t ^''£ ^)' + ' r^;"^} q"(co S ^(co B ^). (8) 
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Comparing with Eq. (|3]), we deduce: 
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This is a new addition theorem for the Legendre function of the second kind. Note that 
Q m _| = Q_ m _i as per Eq. (8.2.2) in @. 

Similarities and contrasts between the pairs of equations, Eqs. ([I]) and (0) and Eqs. ([|) 
and (]^), are worth emphasizing. Of the four variables (r</r>, 6, 9', <fi — </>'), the first pair of 
equations, Eqs. (|T]) and (0), expresses the inverse distance as a series in powers of the first 
variable with coefficients Legendre polynomials of the first kind in 7, which is a composite of 
the other three variables and decomposable in terms of them through the addition theorem 
as in Eq. (||). The second pair of equations, Eqs. (|3]) and (|]), on the other hand, expands 
in Eq. (|j) the inverse distance in terms of the variable <fi — <p', with expansion coefficients 
Legendre functions of the second kind in x, which is a composite of the other three variables 
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FIG. 1: Alternative double summations in (£,m) space (a) first over m at fixed I to form partial 
sums "do as in Eq. (0) (b) first over I at fixed m to form partial sums Lp m as in Eq. (||) 

(r < /r > ,8,8 > ) and decomposable in terms of them through the addition theorem in Eq. (^j). 
For this comparison, it is useful to recast Eq. ([I]) in the more suggestive form, 
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Whereas this expansion has half-integers in the exponents and integer degree Legendre 
polynomials of the first kind as coefficients, Eq. (|]) has integer m's in the exponents and 
half-integer degree Legendre functions of the second kind as coefficients. 

Yet another alternative to Eq. (j3|) follows upon casting the square root in the expression 
for the distance in terms of r, r', and 7 in the form Eq. (El) through the definition 
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This gives the expression 
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now a Fourier expansion in 7 instead of the angle (0 — 0') of Eq. with Q n _i(v) as 
the coefficients. In terms of hyperspherical coordinates, widely used in atomic and nuclear 
study of three (or more) bodies 0, the variable v is esc 2a, a being a "hyp er angle" . In an 
Appendix, we present a number of alternative expressions for the functions Q m _i, some of 
them new, which we have found useful while working with the expansions in Eqs. (|3]) and 

(0)- 

III. TWO-ELECTRON INTEGRALS 

One familiar application of expressions for the inverse distance is to the Coulomb inter- 
action between two charges. We contrast usage of the alternative expansions in Eqs. ([!]) 
and (D) for calculating the electrostatic interaction as it appears in atomic, molecular and 
condensed matter physics. Thus, consider the so-called "direct" part of this interaction 
between two electrons in the 3d 2 configuration, 

v£ = Jdxj d^^m^')^^^)^')- (13) 



The standard treatment |H| uses Eqs. ([!]) and (§), carries out all the angular integrals 
through Racah-Wigner algebra, leaving behind radial "Slater integrals" F k (dd),k = 0,2,4, 
and yielding (for illustrative purposes, all m values have been set equal to zero) 

V£ = F°{dd) + {A/49)F 2 {dd) + (36/441)F 4 (cM), (14) 

where the coefficients are evaluated in terms of Wigner 3j-symbols or are available in tables 
ITJ . The Slater integrals, 

F\dd) = J J r 2 drr' 2 dr\ryr^)RUr)RUr^ (15) 

remain for numerical evaluation. In this example, upon evaluation with hydrogenic radial 
functions, we obtain V = 0.092172 in atomic units. 

The alternative calculation through Eq. (|3|) involves only the m = term and thereby 
the integral 

V e D e=lf J J J ^ 1/2 ^' 1/2 rf J R , ^'Q_i(x)l^(x)| 2 |^(x')| 2 . (16) 

The integrand is a function of z and R variables alone and our numerical evaluation of this 
integral reproduces the value cited above. 



As a second example, we computed an exchange integral for the 3d4f configuration again 
setting, for simplicity, all m equal to zero: 

V£ = J J dxdxV 3 *,(x)^ / (x0r^^V'3 d (x / )V'4/(x). (17) 

The standard method through exchange Slater integrals G and Wigner coefficients gives fll^ 

V e E e = (9/35)G\df) + (16/315)G 3 (d/) + (500/7623)G 5 (d/) (18) 

and, again through hydrogenic radial functions, gives the value V E = 0.0082862. We re- 
produce the same result upon directly computing Eq. (|17D with Eq. (|3]), again involving a 
single four- dimensional integral as in Eq. ( |T6"D with Q_i. 

As the orbital angular momenta involved of the two electrons increase, the number of 
terms in expressions such as Eqs. (|T4| ) and ( [IB] ) also grows, necessitating the computing 
of more Wigner coefficients and Slater integrals. By contrast, only a single term of the 
expansion in Eq. (|3[) and a single integral is necessary in our suggested alternative, the <fi 
integrations setting m = for direct terms and m equal to the difference in the m values 
of the two orbitals for exchange terms. This same selection rule imposed by the <fi inte- 
grations means that even in a calculation with several configurations and the imposition of 
antisymmetrization, such as in a multi-configuration Hartree-Fock scheme, matrix elements 
of |x — x'| _1 between each term in the bra and in the ket gets a contribution from only one 
m value in the expansion in Eq. (|3]). This is a significant economy. 

Although the evaluation of the four-dimensional integrals in Eq. (|T6|) is computationally 
more demanding than the two-dimensional integrals of Eq. ([15]), the preparatory work of 
the Wigner-Racah algebra is avoided. A more efficient alternative to directly evaluating 
the four-dimensional integrals in expressions such as Eq. ( |T6| ) is to use the fact that they 
can be viewed as electrostatic interaction energies between charge densities given by the 
product of wave functions. The charge density p(x') is then used to find the potential at x 
and the resulting potential integrated with the appropriate charge density p(x) to compute 
the integral [1_1|. This same approach applies to the decomposition according to Eq. (|3]) for 



each m-th component of the potential. A two-dimensional Poisson equation in the primed 
variables in Eq. ([TBI) is first solved and then Eqs. (0) and ( |17|) evaluated as two-dimensional 
integrals over the unprimed variables |12|| . As already noted, direct integrals will involve 
only the m = component; exchange ones involve a single m value equal to the difference 
in the azimuthal quantum numbers of the two orbitals. 



IV. SUMMARY 



The inverse distance between two points x' and x is intimately involved in Coulomb and 
gravitational problems. Its expansion in terms of Legendre polynomials Pi of the angle 
between the vector pair or a further double-summation expansion involving the individual 
polar angles of the vectors are well known and widely used in physics and astronomy. We 
have discussed an alternative in terms of cylindrical coordinates, a single summation in terms 
of Legendre functions Q m _i of the second kind in a pair variable \ or double summations 
involving the individual coordinates. These expansions are better suited to problems that 
are decomposable in azimuthal symmetry as shown by applications in |3j and by an illustra- 
tion here for very common electron-electron calculations throughout many-electron physics. 
Further variants are possible for other coordinates such as ring or toroidal, parabolic, bi- 
spherical, cyclidic and spheroidal JL3|, and we hope to return to them in future publications. 



Connections to the theory of Lie groups will also be of interest [14]. 
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APPENDIX A: ALTERNATIVE EXPRESSIONS 

We present in this appendix a number of alternative expressions for the functions Q 
which are useful in calculations using expansions such as Eqs. (|3]) and (O). Setting 6 = 6' 
7r/2 in Eq. @ and using Eq. (8.756.1) of |15|] gives 



n (^- y /r<\M r(l-m + l) n2 2m 
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which can be rewritten as 



Q m x (,) = *e-l"+fr £ 2^" f ^ ) f ^ + 2m 1 ) e^, (A2) 



where we have defined r</r> = e~ v , v = cosh 77. Although the E-th term of these series 
is in different form from what one obtains through the more familiar formula for Q as a 
hypergeometric function |16|], namely, 



Q m _ h (v = cosh??) = ^j^V ^ ^,m+i;m+l ;e - 2 '), (A3) 

their equivalence follows from straightforward algebra. Also, another standard expansion 
for Q in powers of v as in Eq. (8.1.3) of [0, 



is equivalent. However, the results directly in powers of r</r> in Eqs. (|A1|) , ( |A2|) , and (|A3|) 



are more convenient in many applications. Among specific features worth noting in these 
alternative expansions are that only even powers of (r</r>) occur in the sum in Eq. ( |A2|) 
and that for any m, the sum in Eq. ([AT] ) runs over all £ values compatible with it, t > \m\, 
as per their interpretation as angular momentum quantum numbers. 

In the multipole expansion in Eq. ([I]), 7 is an angle formed out of the set (9, 9', — 0') and, 
therefore, cos 7 in the functions P? has range of variation from -1 to 1. On the other hand, 
in the expansions in Eqs. (j^) and (|5[), the arguments v and x °f the Legendre functions of 
the second kind range from 1 to infinity and, therefore, can be written in terms of hyperbolic 
functions as cosh 77 and cosh£, respectively. From Eq. (|j) we have the link between them, 

cosh 77 = cos 9 cos 9' + sin 9 sin 9' cosh £. (A5) 

This disentanglement of v (or if) in terms of a triad is the counterpart of Eq. @ and may 
be used with addition theorems given in the literature such as |5|. || 

00 Vim — n — -) 

Q m ^(coshrj) = £ (-If _ I Q-_x(cosg)P-_x(cosgO e< (A6) 

n=-oo \ 2^ 

An alternative to writing x as cosh£ is to set x = coth^ in Eq. ( |A5| ) . In that case, it follows 
that 



cosh 77 — cosf^ + 9') 



Other variants of the expansions in Eqs. (||) an d (|A6|) follow from identities satisfied by 
the Legendre functions P and Q. In particular, there is an interesting pair of relations 
involving index interchange given in Eqs. (8.2.7) and (8.2.8) of [fj: 



Q£_x (cosh??) = (-1)"* W^-^r(m - n+ I) ^(cothr,), (A8) 



Q^(cothr,) = {-ir f*f V r(m ^ + , } (cosh q) ■ (A9) 
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